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1.  Introduction. 

The  military  and  many  civilian  organizations  require  that  a  suffi¬ 
cient  supply  of  materials  he  in  storage  at  any  one  time  to  meet  certain 
emergency  needs  that  may  arise.  The  fulfillment  of  this  requirement  is 
primarily  a  procurement  problem.  However,  the  procurement  of  an 
adequate  supply  of  materials  alone  does  not  necessarily  guarantee  an 
adequate  supply  of  usable  materials  when  the  need  for  them  arises. 

Their  deterioration  in  storage  may  render  them  unusable.  Despite  all 
precautions,  at  least  a  small  rate  of  deterioration  is  inevitable  for 
all  items  in  storage.  It  is  therefore  essential  that  a  systematic 
program  for  the  surveillance  of  such  stored  items  be  maintained  in  order 
to  insure  that  the  supply  of  usable  items  is  always  adequate. 

The  development  of  such  surveillance  programs  has  received  very 
little  attention  in  the  literature.  Derman  and  Solomon  [1]  studied  the 
quality  of  lots  in  storage  over  time  when  these  lots  are  periodically 
subjected  to  standard  acceptance  sampling  plans.  Ireson  [2]  discussed 
administrative  aspects  of  surveillance  programs.  Savage  [j]  briefly 
examined  when  a  piece  of  equipment  reserved  for  emergency  use  should  be 
withdrawn  from  storage  for  repair;  this  paper  considered  in  general 
terms  the  case  where  the  relevant  costs  are  those  for  regular  repair  and 
for  repair  when  the  emergency  occurs  and  both  these  costs  are  known 


deterministic  functions  of  the  time  elapsed  since  the  last  repair. 

However,  no  comprehensive  investigation  has  yet  been  reported  regarding 
the  appropriate  time  for  corrective  action  to  compensate  for  the  deteri¬ 
orating  quality  of  a  lot  in  storage,  given  the  cost  of  appropriate 
corrective  action  and  the  cost  of  an  imperfect  lot  when  the  emergency 
occurs  as  functions  of  the  number  of  defective  items  in  the  lot.  This 
paper  presents  the  results  of  such  an  investigation. 

The  assumption  made  about  the  lot  in  storage  is  that  the  storage 
life  of  each  item,  i.e.,  the  time  in  storage  until  becoming  defective, 
has  a  common  known  exponential  distribution.  Of  course,  this  condition 
is  not  always  met.  However,  the  exponential  distribution  of  life  has 
been  demonstrated  and  Justified  on  empirical  and  theoretical  grounds  for 
numerous  types  of  items.  Furthermore,  the  particular  exponential  distri¬ 
bution  can  often  be  estimated  by  experimentation  or  from  prior  experience 
with  the  same  product.  Thus,  this  formulation  will  often  become  appli¬ 
cable  at  some  time  after  the  initiation  of  a  continuing  surveillance 
program  for  a  given  product. 

It  is  also  assumed  that  the  time  until  the  emergency  necessitating 
the  use  of  the  lot  has  a  known  exponential  distribution.  This  distri¬ 
bution  should  ordinarily  be  the  appropriate  one  since  it  implies  the 
usual  situation  that  the  emergency  is  a  random  event  in  time.  Examples 
of  such  emergencies  that  might  be  of  interest  sire  the  beginning  of  a 
war,  a  natural  disaster,  the  breakdown  or  failure  of  certain  equipment, 
the  need  to  enter  a  fallout  shelter,  a  rush  order,  the  rejection  of 
regular  material,  and  the  cancellation  of  needed  deliveries  by  a  supplier. 

With  the  information  available  regarding  the  quality  of  the  lot  over 
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time,  it  is  no  longer  necessary  to  sample  the  lot  periodically  to  esti¬ 
mate  its  quality  in  order  to  determine  if  any  corrective  action  is 
advisable.  It  is  now  practical  to  instead  select  the  future  time  when 
the  deteriorating  quality  of  the  lot  will  Justify  corrective  action,  if 
the  emergency  has  not  occurred  yet.  This  corrective  action  might  involve 
inspecting  the  entire  lot  and  replacing  or  repairing  those  items  found 
to  be  deficient.  However,  if  inspection  is  sufficiently  costly,  or  if 
it  involves  destructive  testing,  or  if  the  future  quality  of  the 
non-defectives  in  the  lot  is  mi strusted-de spite  the  assumption  of  an 
exponential  storage  life-  the  proper  action  would  be  to  immediately 
replace  or  repair  the  entire  lot.  The  formulation  presented  here  permits 
both  of  these  alternatives. 

In  the  absence  of  any  surveillance  sampling,  there  remain  only  two 
types  of  costs  that  are  not  independent  of  the  surveillance  program 
selected.  One  type  is  the  cost  of  taking  corrective  action  periodically 
until  the  emergency  occurs.  Its  expected  value  tends  to  decrease  as  the 
time  interval  between  successive  corrective  actions  is  increased.  The 
second  type  is  the  cost  or  the  imputed  cost  incurred  because  of  the 
defective  items  in  the  lot  when  the  emergency  occurs.  The  expected 
value  of  this  type  of  cost  increases  as  the  time  between  corrective 
actions  is  increased.  Therefore,  this  time  needs  to  be  carefully 
selected  so  as  to  obtain  the  proper  balance  between  these  two  types  of 
costs.  To  do  so  requires  that  the  cost  of  corrective  action  and  the 
cost  of  using  an  imperfect  lot,  each  as  a  function  of  the  number  of  the 
defectives  in  the  lot,  be  known. 

The  problem  considered  here  is  the  determination  of  the  time  between 
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corrective  actions  which  minimizes  the  total  expected  cost  from  the 
beginning  of  the  program  until  the  emergency  is  met.  Section  2  presents 
the  formulation  of  the  problem  and  the  operating  characteristics  of  a 
given  surveillance  program.  The  critical  equation  is  derived  for  general 
cost  functions  in  Section  3-  It  is  then  developed  and  simplified  for 
specific  cost  functions  of  interest  in  the  following  two  sections.  These 
results  are  illustrated  with  an  example  in  Section  6.  The  discussion  of 
the  results  in  Section  7  includes  a  suggestion  of  how  they  can  also  be 
applied  in  the  initial  absence  of  information  regarding  the  distribution 
of  storage  life.  The  final  section  summarizes  the  conclusions. 

2.  Formulation  of  the  Problem. 

A  lot  composed  of  N  items  Is  being  held  in  storage.  It  will  be 
used  only  if  certain  emergency  needs  arise,  as  in  the  stockpiling  of 
vital  material.  The  probability  density  function  of  the  time  until  such 
an  emergency  arises,  g(t),  is  given  by 

if  t  <  0 

if  t  >  0  , 

where  0  >  0  and  is  known.  The  quality  of  the  lot  can  be  expected  to 
deteriorate  with  age.  At  periodic  intervals,  corrective  action  is  taken. 
This  corrective  action  consists  either  of  replacing  the  entire  lot  or  of 
inspecting  each  item  in  the  lot  and  either  replacing  or  repairing  those 
items  found  to  be  defective.  The  known  cost  of  taking  corrective  action 
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when  there  are  n  defective  items  in  the  lot  is  R(n).  For  each  item 
in  the  lot,  the  cumulative  density  function  of  the  time  elapsed  since 
the  last  corrective  action  until  the  item  becomes  defective,  F(t),  is 
given  by 


F(t)  =  1  -  e 


-u(t+to) 


> 


-dt 

where  1  -  e  is  the  known  initial  probability  of  being  defective 

and  d  has  a  known  positive  value.  Thus,  at  time  t,  the  defective 
vs.  non-defective  status  of  the  respective  items  can  be  considered  as 
independent  Bernoulli  trials.  Hence,  the  probability  that  the  lot 
contains  n  defective  items  at  time  t,  the  time  elapsed  since  the 
last  corrective  action,  is  given  by 


(I  N\ 
n 


P  (t)  =< 

n  ' 


1  -  e 


-d(t+to)ln 


-d(t+t 


N-n 


,  if  n  =  0,  1,  2,  . . .  ,N 


otherwise . 


The  known  cost  or  imputed  cost  resulting  from  the  presence  of  n 
defective  items  in  the  lot  at  the  time  the  emergency  arises  is  D(n). 
The  objective  is  to  determine  the  value  of  t,  the  elapsed  time 
between  successive  corrective  actions,  which  minimizes  the  total 
expected  cost  from  the  beginning  of  the  program  until  the  emergency 
is  met. 
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3*  Analysis 

let  T  =  (t^,  tg,  ...  )  be  the  sequence  of  positive  numbers  such 
that  'corrective  action  is  taken  at  the  times  t^,  t^  +  tg,  ...  ,  unless 
the  emergency  intervenes.  Let  the  random  variable  E  be  the  time  at 

which  the  emergency  occurs.  Thus,  corrective  action  is  taken  at  the 

m  m  m+1 

times  t^ ,  t  +  t  ,  . . .  ,  Yt  t  ,  where  J  t  <  E  <  ^  t  . 

i=l  i=l  1  "  i=l  1 

The  total  cost  from  the  beginning  of  the  program-! until  the  emergency 
is  met  is 


m 

C(T)  =  ^  R(ni)  +  D(n^)  , 


where  n. ,  n*,  ...  ,n  are  the  number  of  defectives  in  the  lot  at  times 
12'  m 

m 

t^,  t^  +  tg,  ...  ,  £[  ti,  and  n^,  is  the  number  of  defectives  in  the 

lot  at  time  E.  Let  Cq(T)  be  the  cost  from  the  beginning  of  the 
program  until  and  including  the  time  t^.  Let  C^(T)  be  the  cost  of 
the  program  after  the  time  t^.  Therefore, 


C(T)  =  Cq(T)  +  C1(T)  . 

Since 


C  (Tim  =  0)  =  0  , 

CX(T)  =  C^tU  >  0)  P(m  >0)  . 


It  is  assumed  that  the  program  is  begun  immediately  after  corrective 
action  has  been  taken.  Corrective  action  may  be  defined  to  include 
the  initial  placement  of  the  lot  in  storage  if  the  effect  on  the 
quality  of  the  lot  is  indistinguishable. 
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E(C(T) )  =  E{Cq(T)  }  +  E{C1(Tim  >  0)}  P(m  >  0)  . 


I 

Hence, 


Let 

S.  ={t!e{C(T))  =  inf  EfC(X) ),  TeSj  , 

XeS2 

where 

Sg  =  {(t^,  tg,  ...  )  I  tj  >  0,  i  =  1,  2,  ...  } 

i 

It  is  assumed  that  is  not  empty,  i.  e.,  that  there  exists  a 

solution  to  the  problem.  The  objective  is  to  find  a  solution  by 
determining  an  element  of  S  . 

Notice  that 

-at 

P(m  >  0)  =  P(E  >  t  )  =  e  . 

Therefore,  forming  the  expectation  for  Cq(T),  E{C(T))  becomes 

E{C(T)}  =  f  1  (  E  °(J)  P,(t))0e"9t  dt 
Jo  j=o  J 

-9t  n  -et 

+  e  1  E  R(j)P  (t  )  +  e  1  E{C  (Tim  >  0)}  . 
j=o  J 
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Due  to  the  nature  of  the  exponential  distribution,  i.  e., 

P(E  >  t±  +  t! E  >  tx)  =  P(E  >  t)  , 

and  also  because  the  quality  of  the  lot  after  t^  is  independent  of 
its  quality  before  t^,  it  is  clear  that  C^(Tlm  >0),  and  therefore 
EfC^CTlm  >0)},  is  independent  of  t^.  Hence,  it  is  evident  that 
E(C(T)}  is  a  function  of  t^  and  E(C^(Tlm  >  0}},  where  EfC^CTlm  >  OB 
is  a  function  of  tg,  ty  ...  ,  and  that  if  (t^,  tg,  t y  ...  ) 
minimizes  E {C (T) },  i.e.,  if  (t^,  tg,  ty  ...  )eS^,  then 
(t^,  tg,  ty  ...  )  minimizes  E{C^(tI m  >  0) } .  Therefore,  if 

I  t 

(t^,  *2>  )eS]_>  then  if  (t^,  t£,  ty  ...  )  also  minimizes 

EfC^Tlm  >  0)),  (t^,  tg,  ty  ...  )eS^.  It  also  follows  from  the  nature 

of  the  exponential  distribution  and  the  fact  that  the  quality  of  the  lot 

depends  only  on  the  time  elapsed  since  the  last  corrective  action  that 

if  (t^,  tg,  ty  ...  )eS^,  then  (t^,  t^,  tg,  ...  )  minimizes 

E{C. (Tim  >  0)}.  Therefore  if  (t. ,  t t_,  ...  )eSn,  then 
X  X  d  p  X 

(ti,  t^,  tg,  ...  )eS^.  By  repeated  application  of  this  argument,  it 
follows  that  if  (t^  tg,  ty  ...  )eS^,  then  (t^  t^,  t^,  ...  JeS^. 

The  important  conclusion  is  that  one  solution  to  this  problem  is  of  the 
form  that  the  time  between  corrective  actions  is  always  the  same.  This 
is  the  solution  that  will  be  sought.  Thus  denoting  this  constant 
positive  time  between  consecutive  corrective  actions  by  t^,  the 
problem  now  consists  of  determining  the  optimum  value  of  t^. 

It  will  prove  very  useful  to  notice  that,  for  this  newly  defined 
problem,  E{C(T)}  =  E{C^(Tlm  >  0)}.  Since  T  is  now  restricted  to  be 
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of  the  form  (t^  t  ,  t^,  ...  )  where  >  0,  it  will  be  more 
meaningful  to  now  denote  E{C(T))  by  E{;  (t^)J.  The  expression  for  the 
total  expected  cost  can  now  be  written  as 


EfC^)}  = 


+ 


-et 

D(n)  Pn(t))ee 


dt 


R(n)Pn(t1)  +  E{C(t1) ) 


Therefore, 


-0t.) 

i  C 1 

!  N 

1-e  1 

E{C(t  )}  =  1  1 

i  t  D(n)  Pn(t) 

i 

1  A  Jo 

in=o 

-0t 


dt 


-0tn  N 

+  e  1  £  R(n)P  (t  )  , 

n=o 


so  that 

rt.  N  -0t  -0t  N 

I  1  (  C  D(n)Pn(t))0e  dt+e  £ 
o  n=o  n=o 

For  the  particular  P^(t)  function  for  this  problem,  specified  earlier, 
the  conditions  for  Liebnitz's  Rule  are  met  and  the  derivative  exists. 
Setting  the  derivative  equal  to  zero. 


Etc(ti»  fee 

1-e  1 
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=  0  . 


Cancelling  common  terms,  the  critical  equation  becomes 


-tft,l  T  N  ..IK  K 

1-  E  D(n)P„(ti)  +  i  £  E(n)P  (t  )  -  E  H(»)Pa(t1) 

J  |_n=o  1  ^n=o  n=o 


rt  .  N  a.  -9 1,  N 

-  I  1  t  D(n)P  (t)  ■;  9e~dt  dt  +  e  1  j  R(n)P  (t  )  . 
Jo  n=o  n=o 
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Simplifying  slightly. 


N 

E 


n=o 


D(n)Pn(t1) 


+  1  JL 

9  dt. 


N 


£  R(n)Pn(t1) 


D(n)Pn(t)  0e 


-0t 


dt  + 


N 


2  R(n)Pn(t1) 
n=o 


It  is  apparent  that  this  critical  equation  vould  be  exceedingly 
difficult  to  solve  for  t^  for  most  D(n)  and  R(n)  functions. 

But  if  one  were  unable  to  solve  this  equation,  the  analysis  would 
have  little  practical  value.  Fortunately,  greatly  simplified 
equations  quite  amenable  to  solution  can  be  obtained  for  certain  D(n) 
and  R(n)  functions  of  great  practical  interest.  This  is  especially 
true  where  D(n)  and  R(n)  are  linear  functions  of  n.  This  case 
is  studied  and  a  comprehensive  set  of  solutions  graphed  in  the  next 
section.  The  following  section  then  studies  the  case  where  R(n)  is 
a  linear  function  and  D(n)  is  a  certain  quadratic  function. 
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4.  Solution  for  Linear  D(n)  and  R(n)  Functions 


Assume  that  D(n)  =  C^-n.  In  other  words,  the  cost  of  having 

defectives  in  the  lot  when  the  emergency  occurs  is  proportional  to  n, 

the  number  of  those  defectives.  This  would  seem  to  be  the  appropriate 

D(n)  function  for  many  situations.  It  should  also  serve  as  a  suitable 

approximation  for  many  more  situations. 

Assume  that  R(n)  =  K  +  C  *n,  where  K  >  0.  In  other  words,  the 

cost  of  taking  corrective  action-whether  it  be  the  replacement  of  the 

entire  lot,  or  more  likely  the  screening  of  the  lot  and  the  repair  or 

replacement  of  defective  items  -  consists  of  a  fixed  cost  plus  a  certain 

cost  for  each  defective  item  in  the  lot.  This  would  appear  to  be  the 

appropriate  R(n)  function,  or  at  least  a  reasonable  approximation, 

for  practically  every  situation.  The  fixed  cost  could  be  the  red  tape 

cost  of  ordering  or  producing  new  items  plus  the  cost  of  inspecting  the 

lot.  The  cost  for  each  defective  item  could  be  the  cost  associated  only 

with  replacing  or  repairing  that  item.  If  the  lot  is  automatically 

replaced,  the  entire  cost  is  K. 

When  performing  the  summation  operations  for  this  case.  It  is 

crucial  to  notice  that,  for  fixed  t,  the  Pn(t)  probabilities  describe 

a  binomial  distribution  with  parameter  (l-e  and  that  the 

-H(t+t  )  N 

mean  of  this  distribution  is  N(l-e  °  ).  The  D(n)Pn(t)  and 

N  n=0 

R(n)P  (t)  summations  are  analogous  to  finding  this  mean. 

n=o 


Thus, 
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N  N  _  -4 (t+t  )  — M-  ( t+t  )  , 

E  D(n)Pn(t)  =  ^  C  »n(^)  (1-e  °  )n  (e  °  ) 

n=o  n=o 


-u ( t+t  ) 

=  CD*N(l-e  0  )  , 


and 


N  N  -4 (t+t  )  -4 (t+t  ) 

X  R(n)Pn(t)  =  £  (K+CR*n)  (“)  (l-e  °  )n  (e  °  )N"n 

n=o  n=o 


4  (t+t  ) 

=  K  +  C  *N(l-e  °  ) 
K 


Therefore,  after  the  indicated  differentiation  and  integration  operations 
are  performed,  the  critical  equation  becomes 


r  -etii  r 

1-  C  -N  1- 

■  J  ■L' 


^(ti+to) 


i 

+  £cR-N4e 


cd-n 


etT  ft 

n  1  P  O 

1-6  -  e 


l-e 


-t-^u+e) 


K+Cr-N  I l-e 


-u(t+t0)1 
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Combining  like  terms  and  simplifying, 


-vN 


C_*NW 

R 


-nt. 


CR.Ne 


VN 


e 


-t^n+0) 

e 


r  -Ut 

0 

K+C  *N-C  *Ne  0 

R  D 

|i+0 

Simplifying  further  and  dividing  through  by  N, 


-ut. 


r  „  crw  „ 

I  Orf*  .  I 

l  -«*J 

\  CD  +  ~~6~  +  CR 

l  "cd  +  e  +  cd  n+ej 

e 

Ut 


k  „  „  -Kto  l_e_ 

H  CR  "  CDe  u+e 


Finally, 


“M’t, 


r  .  0+8  r  ,  . 
"CD  0  °R 


JLC  -  £  c 

H+0  D  9  UR 


-0t,  i 


Ut  r  -Ut 

_e  °  K+c  -  Ce  0 
"  e  1  N  CR  V 


9 

d+0 
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In  order  to  simplify  notation,  let 


D  = 


j*  1  p  p  _  o 

e 

N  +  CR  '  CD*  I 

uL+ej 

A  = 


-CD  + 


u+e 


if  D  /  0  • 


B  = 


JL  c 

u+e  ud 


K 


if  D  j*  0  . 


Thus,  assuming  D  /  0,  the  critical  equation  can  now  be  written  as 


A  +  Be 


=  1. 


Notice  that 


Therefore,  the  most  convenient  form  of  the  critical  equation  is 


Ae 


-Ht, 


Notice  that  it  has  a  positive  solution  only  if  A  >  1.  While  sin 
explicit  solution  for  t^  is  not  available,  this  critical  equation  can 
be  easily  solved  for  t^  by  numericed.  methods.  This  is  especially  true 
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since 


<  0  ,  for  t  >  0,  A  >  0  . 

Thus,  if  A  >  ~Q~t  the  critical  equation  has  a  unique  positive  solution 
otherwise,  it  has  no  positive  root.  Furthermore,  because  the  left  side 
of  the  equation  strictly  decreases  as  t^  increases,  t^  can  be 
estimated  as  closely  as  desired  by  determining  lower  and  upper  bounds 
to  t^  by  trial  and  error. 

Figures  4.1,  4.2,  and  4.3  present  a  comprehensive  set  of  solutions 
of  the  critical  equation.  For  the  particular  values  of  (A  -  and 
({1+©)  interest,  the  value  of  Ht^  and  therefore  t^  can  be  closely 
estimated  by  interpolating  between  the  given  level  curves. 
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e 


Figure  4.1.  Solutions  of  the  Critical  Equation 


1-1  +  0 


Figure  4.2.  Solutions  of  the  Critical  Equation. 


o 

o 


It  is  of  course  true  that  being  a  positive  root  of  the  critical 
equation  is  a  necessary,  but  not  sufficient,  condition  for  minimizing 
E{C(t^)}.  The  sufficient  condition,  that 


E(c(t1) 3 

I 


>  0 


y 


will  now  be  investigated.  Assume  that  A  >  —q~  ,  D  /  0. 

Tracing  through  the  derivation  of  the  critical  equation,  it  is 
clear  that 


dt 


E{C(t1)  } 


-9t_ 


9e 


1-e 


6tl  2 


-Ht, 


Ae 


1- 


H+0 


-et  i 


-v-t 


-  1  Ne 


Therefore, 


-et. 


dt. 


E(c<ti)1|  -'rVi  |-*‘e 

1-e  1 


-nt  r  -&t  • 
-AHe  1  1-e 


-Ht 


Ne 


+  D 


-et  -9t1  -et  -et  -et 

■(1-e  1r  9  e  _^ge  I‘2(l-e  I)6e  1 

-et.\4 

1-e  x 


\ 

-nt 

Ae  1 

la  "6tl 
~li+0e 


-Ht 


-1  Ne 


Assume  that  t^  is  a  positive  root  of  the  critical  equation.  Hence, 
the  second  term  vanishes.  Additionally, 
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-et. 


Se 


1-e 


■etl2 


-Ht, 


■AHe 


1-e 


i 


-ut 

Ne  °  <  0 


Therefore, 


dt' 


2  (E{C(t1))|  >0 


if  and  only  if  D  <  0.  In  other  words,  if 


K 

N 


+  Cr 


-Ht 


CD6 


Q 

n+e 


> 


* 

then  the  t^  which  is  the  positive  root  of  the  critical  equation 
minimizes  EpC^)].  This  is  further  verified  by  the  fact  that,  for 
this  case, 


E(C(t1) ) 


<  0 


* 

for  0  <  t  <  t^  • 


It  remains  to  investigate  the  other  cases.  If  D  >  0  it  becomes 
quickly  clear  that  this  implies  A  <  which  in  turn  implies 


_d 

dt 


EtCO^)} 


<  0 


for  t^  >  0  , 
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which  Indicates  that  corrective  action  should  never  be  taken,  i.e., 

t1  =  * .  This  is  intuitively  reasonable  since  D  >  0  Implies  that  the 

cost  of  corrective  action  is  high  relative  to  the  cost  of  having 

defectives  in  the  lot  when  the  emergency  occurs. 

l-L+d 

It  is  easily  seen  that  the  case  that  A  <  -g—  can  not  occur  if 
D  <  0.  Therefore,  the  only  remaining  case  is  D  *  0.  For  this  case, 
the  critical  equation  becomes 


+ 


r 

e  r 


■  o  . 


But  this  equation  has  no  positive  root  since 


.  u+e  .  n+e  n 

“cd  “e~  cr  <  T"  D  *  0  ' 


since  K  >  0.  It  is  similarly  verified  that,  for  this  case, 


d 

dt. 


EtCftj)  )j  = 


.c  ♦eaeW*toja!! 

LD  !TCRjHe  /  -St.r2 

1-e 


-»*t. 


,  |jl|  •etil 

1  U+ej e  J 


<  0  for  t^  >  0 

Therefore,  corrective  action  should  never  be  taken. 

The  general,  conclusion  is  that  the  t^  which  minimizes  E{C  (t^) } 
is  the  positive  root  of  the  critical  equation  if  a  positive  root  exists; 
if  the  positive  root  exists,  it  will  be  unique;  if  it  does  not  exist,  t^ 
should  be  infinite,  i.e.,  never  take  corrective  action. 
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5 .  Solution  for  Certain  Quadratic  D(n)  and  Linear  R(n)  Functions 


2 

Assume  that  D(n)  =  CD*n  .  This  should  serve  as  a  suitable 
approximation  for  the  type  of  situation  in  which  the  availability  of 
non-defective  items  in  the  lot  can  partially  compensate  for  the  defective 
items,  so  that  the  imputed  cost  of  having  only  a  few  defectives  in  the 
lot  is  relatively  low  but  this  cost  increases  very  rapidly  as  the  number 
of  defectives  increases.  Assume,  as  before,  that  R(n)  =  K  +  C_*n, 
where  K  >  0. 

Proceeding  as  before,  it  needs  to  be  recalled  that  the  second 
moment  of  the  binomial  distribution  described  by  the  P  (t)  probabil¬ 
ities  is 


(N2  -  N) 


-M-(t+t  )\ 

(l-e  0  ]  .  S 

,  o' 

l-e 

N 

The  V.  D(n)P  (t)  summation  is  analogous  to  finding  this  second  moment. 
u  n 

n=o 

Thus, 


N 


E  D(n)Pn(t) 
n=o 


n=o 


2 

n 


0 


1-e 


-H(t+t  ) 
o 


n 

-H(t+t  ) 

e 

N-n 


-I-1  (t+t  )| 

2 

-^(t+tV 

l-e 

+  N 

1-6  J 

=  NCt 


(N-l> 


-2H (t+t  ) 
o 


+  (l-2N)e 


-H(t+t  ) 


+  N 


Therefore,  proceeding  in  a  straight-forward  fashion  and  simplifying, 
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ft,  n  -et 

I  S  D(n)Pn(t)Se  dt 

Jo  n=o 


=  NCt 


(N-Dee 


-24t 


24 +e 


-4t 


e_ 
4+5 


-t  (4+e)»  |  -t  0 

1-e  A  +  N  1-e 


Hence,  after  considerable  simplification, 


where 


_d 

dt 


-  (Elctt,)))  - 


-et.  -4t 

9e  1  Ne  0 


1-e 


-eMa 


Q  •  Cd(N-1)« 


-4  (t  +2tJ 
o  1 


.  24  “**1 

1  -  2in?e 


_titl  4  “®tl 

+CD(l-2B)e  1-5=8  e 


-lit  ,  6t  |  -lit 

*°»»'  4-  V  °Ee 


St~  +  St~  +  cd(n-i}  e  °(l?4  +  CD(1'2N)(^ 

,  „  O  0  '  '  11 

L  Ne  e 


Thus,  the  critical  equation  is 


Q  =  0  . 
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Assume  that 


is  a  positive  root  of  the  critical  equation.  Thus,  if 


at: 


E{C(t1) ) 


# 

t1=ti 


>  0  , 


* 

then  E{C(t^)}  is  at  least  a  local  minimum. 


at: 


E{C(t1)} 


de 


6t  -tit 
1  Ne  0  dQ 

' — -ut  — 


1-e 


2  dt. 


-  Ne 


■nt 

0 

rl 

1  -et!  2  2  -0*1  2  -»‘i 

[l-e  «  e  *»  e 

-et 

l-e 

-et  4 

l-e  1 

- 

Q 


Hence, 


dt: 


* 

V*! 


-et  -ut 

Qe  1  Ne  °  dQ 

I  -et*l2  dti 

l-e 


* 

tl=tl 


Thus,  if 


dQ 

dt. 


*  >  0  > 

V*! 


then  E(c(t^))  is  a  local,  minimum.  After  considerable  simplification, 
dQ 


dt. 


can  be  written  as 
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d£ 

dt. 


-wt. 


We 


-et. 


l-e 


2Cd(N-1) 


l./'Wi  t  cD- c 


Let  be  the  value  of  t^  such  that 


-W(t  +t' ) 

2Cd(»-1)  [l-e  01  +  0D  -  0R 


w+e 

e 


Then,  clearly,  If  t1  >  0,  then 


M_<0 

dtx  <  U  » 


if  t±  <  tj_  ; 


49_>o 

dt±>  V  ’ 


if  t  >  t  j  . 


It  now  needs  to  be  pointed  out  that  Q  is  continuous  and  that  Q  <  0 
at  t^  =  0.  Therefore,  there  exists  e  >  0  such  that 

(EfC(t1)))<0  if  0  <  t±  <  e  . 

Hence,  since  <0  if  t  >  0  and  t  <  t  '  , 

ut^'  1  IX 

—  (E(C(t1)))<0  if  t±  >  0  and  <  t^  . 
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Therefore,  if  >  0  and 


~  (EfC^)})  =  0  , 

then  ^  >  t|,  and  hence,  >  0.  The  implication  is  that  if  t^ 

is  a  positive  root  of  the  critical  equation,  then  E{'C(t1))  must  be  a 

local  minimum.  Furthermore,  since  >  0  for  t,  >  0  and  t..  >  t', 

dt^  1  1  ]. 

*  * 
t^  must  be  a  unique  positive  root.  Additionally,  t^  can  easily  be 

estimated  as  closely  as  desired.  Merely  proceed  by  trial  and  error 

*  *  dO 

while  remembering  that  tn  >  0  and  t  >  t' ,  and  that  — >  0  for 

x  xx  at1 

* 

t^  >  0  and  t1  >  tj^  .  Thus,  lower  and  upper  bounds  for  t  can 

immediately  be  established  and  then  revised  closer  and  closer  to  the 
* 

true  value  of  t^.  The  one  exception  to  this  procedure  would  arise 
when  there  does  not  exist  even  the  single  positive  root  of  the  critical 
equation.  This  could  happen  if  and  only  if 


(EfCO^)})  <  0 


for  t^  >  0  , 


which  in  turn  occurs  if  and  only  if 


-Ut 


CD(N-l)e 


C  _nt 

-Ht  +  -Ht  +  CD^_1^e  l  2M.+0 


K 


LNe  °  e  ° 


6 


+  cd  (1-2N)  ln5e 


<  o. 
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i.e.,  never  take 


For  this  case,  the  obvious  recommendation  is  t^  ■ 

corrective  action.  The  final  important  implication  of  the  results  is 

* 

that,  if  the  unique  positive  root  t^  exists,  then  since  Q  <  0  if 

0  <  ^  <  t*  and  Q  >  0  if  t^  >  t*,  and  since  (EfCft.^)})  has 

the  same  sign  as  Q  for  t^  >  0,  E(C(t^)}  must  be  a  global  minimum. 
In  summary,  if 


-Ht 


CD(N-l)e 


K 


-Ht 


-Ut 


+  C. 


-Ut 

(N-l)e  0 


LNe 


2n+e 


cd(i-2N)  i5e 


<  o 


*  # 

never  take  corrective  action.  Otherwise,  select  t1  =  t^,  where  t^ 

is  the  necessarily  unique  positive  root  of  the  critical  equation,  Q  =  0. 

* 

This  policy  will  minimize  E{C(t.)}.  Before  solving  for  t^,  determine 
if 


2(N-1)  1 1-e 


-Ht 


+  1 


CD^CR 


0 


If  so,  then  tt-  >  0  for  t  >  0.  If  not,  determine  t',  or  at  least  a 

1  J. 

lower  bound  to  t|,  where 


2Cd(N-1) 


-4*{t  +t«) 

O  1 


> 


*  dQ 

in  which  case  t^  >  t|  and  >  0  for  t^  >  t^.  In  either  case, 

Q  <  0  for  t^  *  max  (0,  t^}.  Then,  letting  max  {0,  t£  3  be  an  initial 
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>  0  for 


lower  bound  for  and  using  the  fact  that  ^ 

I 

t^  >  max  {0,  t^},  proceed  by  trial  and  error  or  by  more  systematic 

* 

numerical  methods,  to  converge  upon  the  true  value  of  t^  to  whatever 
accuracy  is  desired. 

6.  Example 

A  certain  company  in  the  electronics  industry  uses  an  assembly  line 
for  assembling  one  of  its  products.  One  of  the  initial  inputs  of  this 
assembly  line  is  a  certain  kind  of  electronic  tube.  This  tube  arrives 
from  an  earlier  set  of  operations  in  lots  of  100,  usually  at  a  uniform 
rate  of  ten  lots  each  week.  Unfortunately,  due  to  a  number  of  reasons 
such  as  the  breakdown  of  certain  equipment,  an  occasional  lot  is 
prevented  from  being  delivered  when  needed.  It  has  been  observed  that 
the  failure  of  the  lot  to  arrive  on  schedule  is  a  seemingly  random 
event  in  time,  and  that  the  average  time  between  these  events  is  20 
weeks.  The  consequence  of  not  having  any  of  these  tubes  available  is 
that  the  assembly  line  must  shut  down.  To  safeguard  against  this  very 
serious  situation,  a  spare  lot  of  the  tubes  is  kept  in  reserve  at  a 
centralized  storage  station.  Thus,  whenever  the  lot  in  the  regular 
flow  of  production  is  delayed  or  rejected,  the  spare  lot  is  immediately 
removed  from  storage  and  used  in  place  of  the  regular  lot.  The  new  spare 
lot  is  then  produced  during  overtime  work  if  necessary. 

It  has  been  observed  that,  ordinarily,  the  proportion  of  the  tubes 
which  are  produced  defective  is  negligible.  From  previous  experience, 
both  within  and  outside  the  company,  it  is  estimated  that  the  storage 
life  of  this  tube  is  exponentially  distributed  with  a  mean  of  180  weeks. 
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Therefore,  to  oisfesguar-d  against  using  the  spare  lot  when  it  contains  an 
excessive  numWnte  oof  da-  feet  Ives,  It  is  felt  that  the  lot  should  periodi¬ 
cally  be  scree:s««4  and  the  defective  tubes  replaced  with  non-defective 
tubes.  Howeve:j**r,  i  it  is  not  known  Just  how  often  this  corrective  action 
should  be  takeienen so«o  as  to  minimize  the  total  expected  cost. 

It  is  estrJiM  ^ed  t-hat  the  cost  of  inspecting  the  entire  lot  is 
$100.00.  The  i  repl —acenuents  for  the  defective  items  would  be  produced 
and  inspected  o . duri  .  ng  o-vertime  work,  which  would  involve  a  cost  of 
$10.00  per  itene-ei  •  Sine «  there  is  not  enough  time  to  inspect  the  lot 
whan  it  is  callLlM  into  service,  any  defective  tube  present  would  be 
assembled  into  oo the  r  corresponding  item  of  the  product.  The  result  would 
be  that  the  a8S8.«ibI*led  □product  would  not  test  out;  it  would  then  be 
checked,  partial  My  -  disassembled,  the  tube  replaced,  and  reassembled. 

The  total  extram  co««st  i involved  for  each  defective  tube  is  approximately 

$100.00. 

This  data  sals  ■  sufficient  to  permit  the  use  of  the  results  presented 
earlier  to  provoevide  a  solution.  The  cost  functions  clearly  are 
D(n)  =  lOOn  arml  EER(n)  =  100  +  lOn.  Letting  the  unit  of  time  be  the 
week,  H  =  l/l8ll80  send  6  =  l/20.  Since  the  proportion  of  the  tubes 
which  are  produfctad  defective  is  negligible,  tQ  =  0.  Thus,  in  the  ' 
notation  of  the±e  linmear  model  of  Section  4,  A  =  and 

(A  -  jr)  =  ..  .  Ree fenring  to  Figure  4.1,  it  is  seen  that  Ht  should 

lie  between  O..QO5  sand  0.10.  Linearly  interpolating  between  the  two 
curves,  =  =105555*  Therefore,  t^,  the  elapsed  time  between 
consecutive  lnasnfcfcJfclons  of  -the  same  lot  in  storage,  should  be  about 
9.9  weeks. 
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7 •  Discussion  of  Results 


All  of  the  foregoing  results  have  been  based  on  the  assumption  that 
the  particular  exponential  distribution  of  storage  life  is  known. 
Therefore,  the  following  statement  may  at  first  seem  surprising  and 
contradictory.  It  is  that  the  widest  application  of  these  results 
probably  should  be  to  the  situation  where  the  distribution  of  storage 
life  must  be  judged  to  be  unknown,  at  least  at  the  initiation  of  the 
program.  To  understand  why,  it  should  be  recognized  that,  for  this 
situation,  a  program  of  surveillance  sampling  would  or  should  be 
conducted  which  would  have  the  effect  of  gradually  identifying  the 
underlying  distribution  of  storage  life.  Thus,  it  is  entirely  feasible 
that,  by  the  time  corrective  action  should  be  taken,  the  surveillance 
sampling  can  sufficiently  identify  the  underlying  distribution  of 
storage  life  to  signal  the  need  for  this  corrective  action  with  the  use 
of  the  foregoing  results.  In  this  connection,  it  should  be  noticed  that 


when  the  optimum  t  was  being  solved  for,  it  was  really  dt  that  was 

d 

being  found  for  a  particular  ratio  .  Thus,  since  1-e  is 


the  expected  fraction  defective  in  the  lot  at  time  t^,  and  since  dtQ 

-dt 

is  chosen  so  that  1-e  is  the  expected  initial  fraction  defective, 

the  result  of  the  analysis  is  to  determine  the  value  of  expected  fraction 
defective  at  which  corrective  action  should  be  taken.  Hence,  when 
surveillance  sampling  indicates  or  predicts  approximately  this  fraction 
defective  while  estimating  the  corresponding  ratio  g-,  then  corrective 
action  would  appropriately  be  taken  at  the  predicted  time.  In  summary, 
whenever  the  particular  distribution  of  storage  life  is  unknown,  but 
assumed  to  be  exponential,  use  surveillance  sampling  in  effect  to  estimate 
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the  particular  distribution  and  then  apply  the  corresponding  results. 

One  limitation  of  the  results  is  that  the  distribution  of  storage 
life  is  assumed  to  be  exponential.  However,  this  should  at  least  serve 
as  a  useful  approximation  in  a  large  proportion  of  the  cases.  Even  when 
the  distribution  of  storage  life  is  known  to  be  not  nearly  exponential, 
the  results  might  be  enlightening  and  instructive.  It  should  also  be 
recognized  that,  especially  when  corrective  action  will  consist  of 
replacing  the  lot,  only  the  early  part  of  the  distribution  of  storage 
life  may  really  need  to  resemble  the  exponential  distribution. 

Another  limitation  of  the  results  is  that  solutions  are  readily 
obtainable  only  for  a  limited  number  of  D(n)  and  R(n)  functions. 
Fortunately,  among  these  functions  are  the  extremely  important  ones 
studied  in  the  precedi  3  sections.  The  R(n)  function  studied  appears 
to  be  at  least  a  reasonable  approximation  for  most  situations  of  interest . 
The  two  D(n)  functions  studied  appear  to  be  reasonable  approximations 
for  the  two  primary  types  of  situations.  One  of  these  types  is  where 
the  imputed  cost  of  defective  items  in  the  lot  at  the  time  of  the 
emergency  Is  roughly  proportional  to  the  number  of  such  defectives. 

This  would  occur,  for  example,  where  the  results  of  the  use  of  the  items 
in  the  lot  are  relatively  independent.  The  other  primary  type  of 
situation  is  where  the  imputed  cost  of  defectives  when  the  emergency 
occurs  increases  rapidly  as  the  number  of  these  defectives  increases. 

This  would  occur,  for  example,  where  a  few  defectives  in  the  lot  would 
be  easily  compensated  for  by  the  many  non-defectives  but  the  usefulness 
of  the  lot  would  become  seriously  compromised  by  the  presence  of  a 
relatively  large  number  of  defectives.  The  results  derived  for  the 


52 


quadratic  D(n)  function  should  be  useful  for  this  type  of  situation 
either  as  a  preliminary  result  requiring  adjustments  or  as  a  final 
solution.  The  linear  D(n)  function  results  should  similarly  be 
useful  for  the  first  type  of  situation  discussed. 

The  model  does  not  consider  the  possibility  of  the  obsolescence  of 
the  items  in  storage  or  of  the  termination  of  the  time  in  which  the 
relevant  emergency  can  occur.  However,  unless  such  an  occurrence  is 
expected  soon,  it  does  not  appear  that  the  accuracy  of  the  results 
would  be  affected  significantly. 

8.  Conclusions 

The  results  presented  here  should  be  helpful  in  many  cases  in 
determining  when  corrective  action  should  be  taken  because  of  the 
deteriorating  quality  of  a  lot  in  storage .  The  model  does  assume  the 
knowledge  of  considerable  information,  including  the  common  distribution 
of  storage  life  of  the  items  in  the  lot.  However,  it  appears  to  be 
feasible  to  gather  and  refine  this  information  by  surveillance  sampling 
while  simultaneously  applying  the  indicated  results.  Another  important 
limitation  is  the  assumption  that  the  distribution  of  storage  life  is 
exponential.  Nevertheless,  this  assumption  should  be  entirely 
appropriate  in  many  cases  and  a  useful  approximation  in  many  others,  at 
least  for  instructive  preliminary  analysis.  One  of  the  sets  of  cost 
functions  analyzed  and  shown  to  lead  easily  to  a  solution  should  be 
appropriate  for  most  cases  of  interest. 
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